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Abstract 

Contrary to what has been often claimed in the literature, we clarify that the hard photon- 
gluon cross section AaY^^^^{x) in polarized deep inelastic scattering calculated in the gauge- invariant 
factorization scheme does not involve any soft contributions and hence it is genuinely hard. We 
show that the polarized proton structure function up to the next-to-leading order of Og is 

independent of the factorization convention, e.g., the gauge-invariant or chiral-invariant scheme, 
chosen in defining Aa^^^^{x) and the quark spin density. Thereby, it is not pertinent to keep 
disputing which factorization prescription is correct or superior. The hard-gluonic contribution to 
T^, the first moment of Qiix), is purely factorization dependent. Nevertheless, we stress that even 
though hard gluons do not contribute to in the gauge-invariant scheme, the gluon spin component 
in a proton, which is factorization independent, should be large enough to perturbatively generate 
a negative sea polarization via the axial anomaly. We briefiy comment on how to study the 
evolution of parton spin distributions to the next-to-leading order of QCD in the chiral-invariant 
factorization scheme. 
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1. Recently we have analyzed the EMC [1] and SMC [2] data of the polarized proton 
structure function gi{x) to extract the polarized parton distributions to the next-to-leading 
order (NLO) of QCD by assuming that = (Q^) = 10 GeV^ for each x bin of the gi{x) 
data [3]. Our analysis is performed in two extreme factorization schemes: gauge-invariant 
and chiral-invariant ones. In the former factorization, hard gluons do not contribute to Ff, 
the first moment of gi{x), and the quark net helicity is dependent but has a local gauge- 
invariant operator expression, whereas gluons do make contributions to Ff and the chiral- 
invariant quark spin does not evolve in the latter scheme. However, we have stressed that 
physics is independent of the choice of the factorization prescription; the size of the hard 
gluonic contribution to Ff is purely a matter of convention, as first realized and strongly 
advocated by Bodwin and Qiu [4] sometime ago. 

Because of the availability of the two-loop polarized splitting functions APj;j\x) very 
recently [5], it becomes possible to embark on a full NLO analysis of the experimental data 
of polarized structure functions by taking into account the measured x dependence of at 
each X bin. Two of such analyses are now available. A NLO QCD analysis is carried out 
in [6] in the conventional MS scheme, a gauge-invariant factorization, within the framework 
of the radiative parton model. While a sizeable and negative sea polarization is required 
to describe all presently available data, the gluon spin density is found to be rather weakly 
constrained by the data. On the contrary, a NLO fit presented in [7] in the chiral-invariant 
scheme shows that the gluon contribution is large and positive. Furthermore, the authors 
of [7] have criticized the gauge-invariant scheme that it is pathological and inappropriate 
because soft contributions are partly included in the hard coefficient function rather than 
being factorized into parton spin densities. More precisely, the hard gluon-photon cross 
section has the expression [4] 
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in the gauge- invariant factorization scheme with //fact a factorization scale. An objection to 
this scheme has been that the last term in Eq.(l) proportional to 2(1 — x) appears to arise 
from the soft region ~ << Aq^j-,, and hence it should be absorbed into the polarized 
quark distribution [7-9]. 

Although the issue of whether or not gluons contribute to F^ was resolved five years ago 
that it depends on the factorization convention chosen in defining the quark spin density and 
the hard cross section for the photon-gluon scattering [4,10], [| the fact that the interpretation 
of Ff is still under dispute even today and that many recent articles and reviews are still 
biased towards or against one of the two popular implications of the measured gi{x), namely 
sea quark or anomalous gluon interpretation, demands a further clarification on this issue. 
In the present paper, we will point out that, irrespective of the soft cutoff, the 2(1 — x) 

^See also a very nice roundtable summary at the Polarized Collider Workshop (University Park, PA, 
1990) on the theoretical interpretation of the measured gi{x) [11]. 
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term in Eq.(l) actually arises from the axial anomaly, i.e., from the region where k'j_ ~ fJ''{acf 
Consequently, many criticisms to the gauge-invariant factorization scheme are in vain. We 
will show that even though hard gluons do not contribute to in the gauge-invariant 
scheme, the bulk of negative sea polarization must be perturbatively generated by gluons 
via the axial anomaly. We would like to emphasize that none of the results presented in 
this paper (it reads like a status review) are new and they are scattered in the literature. 
Nevertheless, we believe that the present paper is useful for clarifing several confusing issues 
and for making a unified and consistent understanding of gi{x). 



2. As far as the first moment of gi{x) is concerned, the parton-model and OPE ap- 
proaches are equivalent. However, in order to consider QCD corrections to gi{x) itself, we 
shall consider the parton model. To NLO, the polarized proton structure function is given 

by 
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'ii^,Q') = ^Y.^^{[Aq,ix,Q') + ^^Af,{x)^Aq,ix,Q')] 
+Aal^^{x,Q')^AG{x,Q')}, 



(2) 



where is the number of active quark flavors, Aq{x) 
AG{x) — G^x) — G^{x), and denotes the convolution 



q^{x) + q^{x) — q^{x) — q^{x), 



fix) g{x) = 



dy 
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X 



9{y)- (3) 



The Afq[x) term in Eq.(2) depends on the regularization scheme chosen, but its first moment 
is scheme independent at least to NLO: Jq Afq{x)dx — —2. Since the polarized photon- 



gluon cross section Aa'^'-^{x) has infrared and collinear singularities at 



and 



= 0, where m is the quark mass, p"^ is the gluon momentum squared and k_i_ is the 
quark's transverse momentum perpendicular to the virtual photon direction, it is necessary 
to introduce a soft cutoff. Depending on the infrared regulators, one obtains 
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(5) 

(6) 



for the momentum regulator (p^ ^ 0) [12], the mass regulator (m^ ^ 0) [13], and the modified 
dimensional regulator (//^ ^ 0) [14], respectively. Note that the term {2x — 1) in Eqs.(4- 
6) is nothing but the spin splitting function 2AP^g(x) and that the term proportional to 
2(1 — x) in (5) and (6) is an effect of chiral symmetry breaking: It arises from the region 
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where fc^ ~ in the mass-re ffulat or scheme, and from /c^ ~ in the n ^ 4 (n > 4) 

MS / \ / 

dimensions in the dimensional regularization scheme. 

The cross sections given in (4-6) are, however, not the desirable perturbative QCD results 
since they are sensitive to the choice of the regulator. Although the ln(Q^/ — p^) and 
InlQ"^ /m^) terms, which depend logarithmically on the soft cutoff, make no contributions to 
the first moment of gi{x), it is important to have a reliable perturbative QCD calculation 
for Aa'^'-'{x) since we are intertested in QCD corrections to gi{x). To do this, we need to 
introduce a factorization scale //fact; so that 

Aa^«(x, Q') = Aal^x, Q\ /^Lt) + Aaj5,(x, /^Lt) (7) 

and the polarized photon-proton cross section is decomposed into 

i 

That is, the hard piece of Ao'^'^ix) contributes to g\[x\ while the soft part is factorized 
into the nonpcrturbativc quark spin densities Aqiix). Since Ao'^'^ix) is a physical quantity, 
a different factorization scheme amounts to a different way of shifting the contributions 
between Aa^^^{x) and Aq{x). An obvious partition of Au^'~^{x) is that the region where 
A*fact contributes to the hard cross section, whereas the soft part receives contributions 
from /c^. ^ A*fact ^^'i hence can be interpreted as the quark and antiquark spin densities 
in a gluon, i.e., Aal^^^{x^[}\^^ = Aq^{x, iif^^). Physically, the quark and antiquark jets 
produced in deep inelastic scattering with kj_ <^ fi^^^^ are not hard enough to satisfy the jet 
criterion and thus should be considered as a part of one-jet cross section [12]. The choice 
of the "ultraviolet" cutoff for soft contributions specifies the factorization convention. In 
the present paper we only focus on two extremes: the chiral-invariant scheme in which 
the ultraviolet regulator respects chiral symmetry, and the gauge-invariant scheme in which 
gauge symmetry is respected but chiral symmetry is broken by the cutoff. 

For a massless quark, one will expect that, based on helicity conservation or chiral symme- 
try, the quark spin Aq = Aq(x)dx is independent and that there is no sea polarization 
perturbatively induced by hard gluons. One way of calculating Aa]^{x) is to make a direct 
cutoff on the k± integration so that its integral expression is exactly the same as Aa'^'^{x, Q"^) 
except that k^^ is integrated over from to (5^(1 — x)/4x for the latter, but from to //fg^^ 
for the former. For fi^^^ » m^, — p^, the results are [4,15] (for a complete expression of 
Aa]^^{x) to NLO using mass or momentum regulator, see [16]) 



^'^iSti^^fJ'Ucthi = Aqg{x,fJ,U = ^ 2" 
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' (9) 



for various soft cutoffs, and the subscript CI indicates that we are working in a chiral- 
invariant factorization scheme. Note that, as stressed in [9], the soft cross sections or quark 
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spin densities in a helicity + gluon given by (9) do not make sense in QCD as they are derived 
using perturbation theory in a region where it does not apply. Nevertheless, it is instructive 
to see that Ag^j = Jq Aq^j^{x)dx vanishes when = or — >> m^, as expected. Hence, a 
sea polarization for massless quarks, if any, is produced nonperturbatively. Now it does make 
sense in QCD to subtract AcrJ^^ from Aa'^'-' [see Eqs.(4-6)] to obtain a reliable perturbative 
QCD result for Aa^^: 

Aa£d(^, Q', /^Lt)ci = ^i2x - 1) fin 4^ + In - l) , (10) 

which is independent of the infrared regulators as long as fif^^^ » /i^ . m?. —p^. It is also 
clear that the soft 2(1 — x) term in (5) and (6) drops out in Aa'^^^^{x). Therefore, 

A^hfrd(<5^/^fact)cI = / c/xAa;[f^d(x,Q^/x^^,Jci = (11) 

Jo /vr 

Since gauge invariance and helicity conservation in the quark-gluon vertex are not broken in 
the chiral-invariant factorization scheme, it is evident that Agci does not evolve, consistent 
with the naive intuition, and it corresponds to a nucleon matrix element of a gauge-invariant 
but nonlocal operator (see e.g., Eq.(4.7) of [4]). 

Nevertheless, we do have freedom to redefine Aa^^^{x) and Ag(x) in accord with Eq.(8). 
The fact that there is no gluonic operator at the twist-2, spin-1 level in the approach of 
OPE [17] indicates that there must exist a factorization scheme in which hard gluons make 
no contribution to the first moment of gi{x) and that Ag can be expressed as a nucleon 
matrix element of a local gauge-invariant operator. In this scheme, gluons can induce a sea 
polarization even for massless quarks. This can be implemented as follows. For < /ifact? 
the box diagram for photon-gluon scattering is reduced under the coUinear approximation 
for the quark-antiquark pair created by the gluon to a triangle diagram with the light-cone 
cut vertex 7'^7^ combined with a trivial photon-quark scattering [4]. As a result, Ag'^(x) can 
be also obtained by calculating the triangle diagram with the constraint < fif^^t- order 
to have a non-vanishing Ag"-^ even for massless quarks, evidently we need to integrate over 
kj_ from to oo to achieve the axial anomaly and hence chiral-symmetry breaking [12], and 
then identify the ultraviolet cutoff with /ifact- We see that the desirable ultraviolet regulator 
must be gauge-invariant but chiral-variant owing to the presence of the QCD anomaly in 
the triangle diagram. Obviously, the dimensional and Pauli-Villars regularizations, which 
respect the axial anomaly, are suitable for our purpose. It is found [18] | 

Aggi(x,/iLt) - Aggi(x,/iLt) = - — (1 - x) (12) 

TT 

for /ifact >> ^Ms'^"^' ^^^^^ subscript GI designates the gauge invariant scheme. 
It should be accentuated that the last term in (12) originating from the axial anomaly [12] 

^For a complete expression of AgQj(x), see [16]. Note that the r.h.s. of Eq.(9) in [16] is too large by a 
factor of 2. 
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comes from the region ~ Aij^j^^. As noted in passing, the quark spin distribution in a 
gluon cannot be rehab ly calculated by perturbative QCD; however, the difference between 
Aq^j{x) and Aq^^{x) is trustworthy in QCD. It is interesting to see from Eq.(12) that 

A?gi(/^Lt) = -^^^ for massless quarks, (13) 

that is, the sea-quark polarization perturbatively generated by helicity + hard gluons via 
the anomaly is negative ! It follows that the hard cross section has the form 



1-x 



2^ 




2 , ... i +2(l-a;) 

/^fact ^ / 



(14) 



and its first moment vanishes: 

'' dxAa£,i(x,Q^)Gi = (15) 



/ 

Jo 
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in the gauge-invariant factorization scheme. Bodwin and Qiu [4] have shown generally that 
the gluonic contribution to vanishes so long as the ultraviolet regulator for the spin- 
dependent quark distributions respects gauge invariance and the analytic structure of the 
unregulated distributions. From Eqs.(9) and (12) we also see that, contrary to what has 
been often claimed in the literature [7-9], though the 2(1 — x) term in (5) and (6) drops out 
in Aa'l^j.^{x)ci because it arises from the soft region A;^. ~ rn^.y? . it emerges again in the 
gauge invariant scheme due to the axial anomaly and this time reappears in the hard region 
^1 ~ A*facf Therefore, the hard gluonic coefficient function is genuinely hard ! 

We wish to stress that the quark spin density Aq'^iyX) measures the polarized sea quark 
distribution in a helicity + gluon rather than in a polarized proton. Consequently, Aq'-'{x) 
must convolute with AG{x) in order to be identified as the sea quark spin distribution in a 
proton: 

A?f (^, /^Lt) - Agf (x, /xLt) = -^(1 - x) ® AG{x, /.Lt)- (16) 

TT 

This relation can be also derived from Eqs.(8) and (14) by noting that Aa"'''{x) = e'^S{x — 
Q'^/2p ■ q) [19]. Now it becomes clear from Eqs.(2), (14) and (16) that despite of the 
factorization- scheme dependence of Aq{x) and Aal[^^.^{x) , the physical quantity g[{x) remains 
to he scheme independent up to NLO, as it should be. The choice of factorization is thus a 
matter of convention. Since the valence quark spin distribution Aq^^x) = Aq{x) — Aqs{x) is 
factorization independent, it follows that [19] 

Aqoiix, /x^g^J - Mci{x, n'j^) = ^(1 - x) (g) AG{x, /x^^^J, (17) 

TT 



which leads to 



AqaiiQ') - Aqo^m = -^^AG(g2), (18) 
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where we have set ^f^^^ = Q^. Eqs.(14) and (17) provide the necessary relations between the 
gauge-invariant and chiral- invariant factorization schemes. For a given AG{x, Q"^), the quark 
spin densities in these two different schemes are related to each other via (17). It should 
be remarked that in spite of a vanishing gluonic contribution to in the gauge-invariant 
scheme, it never means that AG vanishes in a polarized proton. Quite opposite to the 
naive expectation, if there is no sea polarization in the chiral-invariant scheme, then the 
size of the gluon spin component in a proton must numerically obey the relation AG((5^) = 
— (27r/as((5^))AgP(Q^) in order to perturbatively generate a negative sea-quark polarization 
AgP(Q^) via the QCD anomaly. In other words, even gluons do not contribute to Fj*, the 
gluon spin can be as large as 2.5 for Aqf^ = —0.10 at = 10 GeV^ provided that Ag^ = 0. 
Recall that the gluon polarization induced from quark's bremsstrahlung is positive (see the 
first moment of APcqix) in Eq.(25) below). 

Phenomenologically, one has to specify the factorization scale /Xfact in order to extract 
the quark and gluon spin distributions from polarized DIS data as the hard photon-gluon 
cross section is dependent of /ifact- In practice, fj,f^^ can be fixed to be (Q^), the average 
of the gi{x) data set. Of course, one should take into consideration the logarithmic term 
ln{Q'^ / fif^^^) in A(Th^^(a;) to fully account for the measured x dependence of at each x 
bin. 

One may ask how to accommodate the aforementioned two different factorization schemes 
in the approach of OPE ? An examination of this issue also provides a clear picture on the 
differences between Agci and Agci- Since the flavor-singlet axial- vector current has an 
anomalous dimension first appearing at the two-loop level [25], the quark spin 

Agci = {p\qii,i5q\p)s^, (19) 

with a spin 4- vector, is gauge-invariant but dependent. The evaluation of the nucleon 
matrix element of involves connected and disconnected insertions (see e.g., [23]). The 
connected and disconnected insertions are related to valence quark and vacuum (i.e., sea 
quark) polarizations, respectively, and are separately gauge invariant. Thus we can make 
the identification: 

{p\J^\p) = (p| Jrb)con + (p| Jrb)dis = T.i^q^' + (20) 

Interestingly, lattice QCD calculations of Ag^^ and Ag^^ became available very recently 
[23,24]. It is found that Aus = Adg = As = — 0.12±0.01 from the disconnected insertion [23]. 

■^A sea-qaurk interpretation of with As ~ —0.10 ± 0.03 at — lOGeV^ [20] has been criticized on 
the ground that a bound |As| < 0.0521q'q52 [21] can be derived based on the information of the behavior 
of s{x) measured in deep inelastic neutrino experiments and on the positivity constraint. First of all, this 
argument is quite controversial [22]. Second, one can always find a polarized strange quark distribution with 
As ^ —0.10 which satisfies positivity and experimental constraints [3]. Third, a sea polarization of order 
—0.11 is also found by lattice calculations [23,24]. 
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This empirical SU(3)-flavor symmetry within errors imphes that the sea-quark polarization in 
the gauge-invariant scheme is indeed predominately generated by the axial anomaly. Recall 
that sea contributions in the unpolarized case are far from being SU(3) symmetric: d > u > s. 
In order to connect to the chiral-invariant scheme, one can write 

= -Kf' + K^'= + K^", (21) 

with = (a,n//27r)e'^^''Me(apA^ - \gfabcA''pAl) and 60123 = 1- Though neither nor 
K*^ is gauge invariant, their matrix elements can be identified with quark and gluon spin 
components in the light-front gauge — [12] (it is not necessary to specify the coordinate) 

s+A^ci = (p|J5+Ip)a+=o, s+AG = (^\{E X A)+b)A+=o = (p|i^+|p)A+=o- (22) 

The quark spin Agci docs not evolve as the current Jg* is conserved in the chiral limit. Of 
course, both AG'(x) and Agci(''^) (not just their first moments!) also can be recast as matrix 
elements of a gauge- invariant but nonlocal operator [4,15], as noted in passing. Applying 
(22) to the axial-current matrix element leads to 

ip\J^\p) = {p\J^\p)con + {P\j^\p)dis + {p\K>Ui. 

E(Agf + AgP-^AG).^ (23) 

where use of Aq^^ = Aq^^ has been made. It is clear that (20) and (23) are equivalent owing 
to the relation (18). 

3. The dependence of parton spin densities are governed by the Altarelli-Parisi 
equations: 



f^Aq^s{x,t) = ^APf/(x) ® Aq^s{x,t), 

d ( Aqs{x, t) \ _ as{t) ( APUx) 2nfAPgGix) \ / Agg(x, t) 



dt\AG{x,t)J 271 \APG,ix) APgg{x) J \AG{x,t) 



(24) 



where AgNs(2^) = Aqi{x) — Aqj{x), Aqs{x) = J2iAqi{x) and t = ln{Q'^/Ql). The complete 
polarized splitting functions up to NLO, AP{x) = AP^^^x) + ^AP^^\x), have been calcu- 
lated in the MS scheme recently [5]. Hence, the NLO evolution of spin parton distributions 
in the gauge-invariant factorization scheme is completely determined. Explicitly, the AP 
equation for the first moment of fiavor-singlet parton spin densities reads [5] 

I /ASGi(t)\ ^ ( t(-2%) \ MSGi(t)\ 

dt[AG{t)) 271 [2 + l^{25 - luf) f + f f ; 1, AG(t) j' ^ ^ 

where ASgi(0 = Jq dxAqs{x, t), /3o = H - and /3i = 102 - f n/. It is clear that AEqi 
to NLO is Q'^ dependent. 
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One may choose to work in the chiral-invariant factorization scheme, so that AEci does 
not evolve with Q'^ . This requires that 

Aa£,(g^)ci = -g, Tjf'^^f Ai^f(x)dx = 0. (26) 

Using the results obtained in the MS scheme, one may introduce a modification on NLO 
anomalous dimensions and hard coefficient functions to transfer from the Gl scheme to the 
CI prescription [26] . However, this transformation cannot be unique since it is only subject 
to the constraints (26). Indeed, three different scheme changes have been constructed in [7]. 
As a consequence, the NLO evolution of polarized parton distributions in the CI scheme 
obtained in this manner [7] is ambiguious and not quite trustworthy as it depends on the 
scheme of transformation. 

We can avoid the aforementioned ambiguities and complications by working in the context 
of the gauge-invariant scheme where NLO polarized splitting functions are known. Once 

the cvohition of the spin parton distributions A^giI^;, Q^) and AG(x, Q^) is determined 
from the AP equations, Agci(a;,Q^) in the chiral-invariant prescription is simply related to 
AgGi(a;, Q^) and AG'(a;, Q^) by Eq.(17). One can check from Eq.(25) that rf(ASci(t))/rft = 
to NLO, as it should be, where AEci(t) = AEGi(t) + {n fUs/ 271)1^0 {t). In the absence 
of direct calculations of NLO polarized splitting functions in the chiral-invariant scheme, 
we believe that this is the right approach for studying the Q"^ evolution of parton spin 
distributions. 

4. To conclude, contrary to what has been often claimed in the literature, we have clar- 
ified that the 2(1 — x) term in A(T^^j(a;) in the gauge- invariant factorization scheme arises 
from the region ~ ii^^^^ and hence is a genuinely hard contribution. We have shown explic- 
itly that the physical quantity gi{x) is independent of the choice of factorization convention. 
The sea-quark interpretation of in the gauge-invariant scheme or the anomalous gluon 
interpretation in the chiral-invariant scheme is purely a matter of factorization convention 
chosen in defining Ag(a;) and ^a'^^^{x). We have emphasized that even though hard gluons 
do not contribute to V\ in the gauge-invariant scheme, the gluon spin component in a proton 
should be large enough to perturbatively generate a negative sea polarization via the axial 
anomaly, recalling that AG(x, Q^) is factorization independent. 

As far as gi{x) is concerned, both GI and CI factorization schemes are on the same footing. 
Thus it does not make sense to keep disputing which factorization prescription is correct or 
superior. Of course, once a set of Aq'Gi(x), AG(x), A(Th^(j(x)Gi or of Ag'ci(x), AG(a;), /S.a^^^{x)ci 
is chosen, one has to stick to the same scheme in all processes. 

In practice, it appears that the use of ^qci^x) is more convenient than ^qci[x). First of 
all, A^Gi corresponds to a nucleon matrix element of a local and gauge-invariant operator, 
and its calculation in lattice QCD became available recently. For Aq'ci, one has to compute 
the matrix element of in the light-front gauge, which will involve much more lattice con- 
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figurations. Second, NLO polarized splitting functions liave been determined very recently 
in the gauge-invariant scheme, and it is straightforward to study the evolution of AqQi{x, Q^) 
through AP evolution equations. We have argued that a NLO analysis of polarized DIS data 
should be first carried out in the gauge-invariant factorization scheme and then related to 
the chiral-invariant prescription, if desired, via Eq.(17). 
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